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The Green-Kubo relation for two models of granular gases is discussed. In the Maxwell model 
in any dimension, the effective temperature obtained from the Green-Kubo relation is shown to be 
frequency independent, and equal to the average kinetic energy, known as the granular temperature. 
In the second model analyzed, a mean-field granular gas, the collision rate of a particle is taken to 
be proportional to its velocity. The Green-Kubo relation in the high frequency limit is calculated 
for this model, and the effective temperature in this limit is shown to be equal to the granular 
temperature. This result, taken together with previous results, showing a difference between the 
effective temperature at zero frequency (the Einstein relation) and the granular temperature, shows 
that the Green-Kubo relation for granular gases is violated. 
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I. INTRODUCTION 

Fluctuation-dissipation (FD) relations in non- 
equilibrium systems have been the subject of consider- 
able attention recently [J, [3, i, 0, i, [1, 0, Hi- In this 
context, granular matter serves as a perfect example 
of a non-equilibrium system. Such systems consist of 
macroscopic grains undergoing inelastic collisions. In a 
cooling state, no energy is injected into the system, and 
the total kinetic energy decreases steadily, approaching a 
trivial steady-state where all grains are at rest. If energy 
is constantly supplied to the system a non- equilibrium 
steady-state (NESS) is reached. This paper focuses on 
the frequency-dependent FD relations in this NESS. 

If the density of grains is low, the system is known as 
a granular gas [l^. In this case, spatial correlations 
are small, and mean-field models, in which particles col- 
lide without regard to their positions, may provide good 
approximations of the system at hand. In such models, 
interactions are governed by a stochastic collision rule, 
with some interaction rate. Two such models are studied 
in this paper, differing in their collision rules. The sim- 
pler of the two is the inelastic Maxwell model [ll|i [l3 1 in 
which the collision rate is taken to be constant for each 
pair of particles in the system. In the second, more real- 
istic model, which we refer to as the mean-field granular 
gas (MFGG), the collision rate of each pair of particles 
is proportional to the particles' relative velocity. 

The FD theorem [isl . [T3 | connects the linear response 
of an equilibrium system to its correlations. Stated in the 
frequency (w) domain, it relates the response function to 
a weak harmonic force to the Fourier-Laplace transform 
of the auto-correlation function of the observable con- 
jugate to that force. The theorem states that, remark- 
ably, the two functions have exactly the same frequency- 
dependence, and moreover, that their ratio is equal to 
the temperature of the system. Since the full frequency- 
dependent functions are not always easy to compute, 
their a; ^ limit is often calculated; the connection be- 
tween fluctuation and response in this limit is known as 



the Einstein relation. 

In a NESS, the response to a weak perturbation from 
the steady-state and the corresponding auto-correlation 
may be similarly defined. Since granular gases are out of 
equilibrium, however, the FD relation is not expected to 
generally hold, and the two functions may have different 
frequency dependencies. If we nevertheless define a FD 
temperature Tpo as the ratio of the correlation and re- 
sponse functions [l|, [l5| , we would expect this quantity 
to depend on frequency. 

This having been said, there is evidence that the FD 
relation may nonetheless closely hold for driven granular 
gases. Molecular dynamics simulations of 2D granular 
gases indicate that TpD is approximately frequency- 
independent, and moreover, is close to the average kinetic 
energy per particle, known as the granular temperature 
Tg- Einstein relations were calculated for the MFGG 
model, yielding a value for TpD different from Tq [ifj . 
Numerically, however, this deviation was smalL and con- 
sistent with direct Monte-Carlo simulations [Iq that ob- 
served no violation of the FD relation (dependence of 
Tp£) on Lo) to within their accuracy. Experiments on air 
fluidized particles have shown that at low densities the 
effective temperature measured from the Einstein rela- 
tions is close to the granular temperature [l^. For the 
ID Maxwell model with a Gaussian thermostat, the Ein- 
stein relation was calculated [l^j and yielded TpD — Tq- 
Additional transport coefficients, calculated in pT], ver- 
ify the Einstein relation in the Maxwell model. Subse- 
quently, we showed that the full frequency-dependent FD 
relation holds exactly for the Maxwell model, driven by 
either a Gaussian or a stochastic thermostat This 
result was recently shown to also hold in 2D with a Gaus- 
sian thermostat [23| . Fluctuation-dissipation relations in 
denser systems, in which correlations play a role, are dis- 
cussed e.g. in [H, m m, m, [131. Time dependent FD 
relations in other systems are studied in d, |2^ . 

I n 1221 we also discussed a simple model of a granular 
gas [23], for which the FD relations were calculated ana- 
lytically. It was shown that in measurements where cor- 
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relations between different degrees of freedom appear, the 
FD relation is violated. In [23| the role of correlations was 
studied in simulations of granular fluids, and was shown 
to affect the violation of the FD relation. However, FD 
violation has not yet been observed in simulations of di- 
lute granular gases, due to their limited accuracy. In this 
paper we provide theoretical reasoning for why FD rela- 
tions should be only slightly violated even in the dilute 
limit. 

We first show that the frequency dependent FD rela- 
tion holds in the Maxwell model in any dimension and 
with both customarily used thermostats. We then show 
that the frequency-dependent FD relation is violated in 
the MFGG, where collision rates depend on velocity. 
Quantitatively, the violation is expected to be small, in 
accordance with previous simulation results. 

II. FLUCTUATION-DISSIPATION RELATION 

The Green-Kubo (GK) formulation fl3| employs the 
frequency-dependent autocorrelation function, defined as 
the Fourier-Laplace transform of the velocity autocorre- 
lation function, 

D{CJ)^- (v(0).v(t)>e— (1) 

where d is the system dimension. Vectorial quantities are 
written in boldface font (e.g. v for the velocity). Note 
that D{lo) reduces to the diffusion coefficient in the zero 
frequency limit. 

The frequency dependent mobility is calculated by ap- 
plying a weak periodic force ^(t) = ^oe*'^* to a single 
particle i and measuring the amplitude A {uj) of its re- 
sulting velocity, averaged over the steady-state: (v^) ~ 
A(tj)e"^*. Then and A are directed along the same 
direction, and the mobility ^{lo) is defined so that 

= A (w) . (2) 

Equivalently, if some general time-dependent force $,{t) is 
applied to particle i, and (v^ {t)) is its resulting average 
velocity, /i {lo) may be defined so that 

IHmH = (v.H), (3) 

where ^ (cj) and (v^ (w)) are the Fourier transforms of 
^(i) and (v^ [t)) respectively. Note that both D(ijo) and 
//(w) are complex. 

In thermodynamic equilibrium the FD theorem, or 
Kubo formula, guaranties that D{uj) = Tfi{uj), where the 
temperature is measured in units of energy, and Boltz- 
mann's constant is set to one. The Einstein relation is 
obtained by taking the w — > limit, thereby connecting 
the mobility with respect to a constant force and the long 
time diffusion coefficient. 

For a NESS one may use the fluctuation D{uj) and 
the response fi{uj) to define an effective FD temperature 
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Tpd{uj) scales fluctuations in the system and in general 
depends on the measurement frequency to. In this con- 
text, the FD theorem states that in equilibrium TpD {us) 
is independent of uj and is equal to the system's temper- 
ature. 



III. INELASTIC MAXWELL MODEL 

In this section we show that in the Maxwell model, 
where the collision rate F is velocity independent, the 
frequency dependent GK relations hold exactly with 
TpD = Tq, where Tq, known as the "granular tempera- 
ture" , is defined as Tq = (v^) /d (taking particles of unit 
mass). In [20j| the Einstein relations for the ID Maxwell 
model with a Gaussian thermostat were shown to give 
Tfd (0) = Ta (the thermostats are defined below). In 
[2^ we showed that for both Gaussian and stochastic 
thermostats, in the ID Maxwell model Tpu (w) — Tq, 
independent of frequency. For clarity of presentation, 
before turning to higher dimensions and to the more re- 
alistic MFGG model, wc first present the derivation of 
this result in more detail. 

In the ID inelastic Maxwell model, each particle is de- 
scribed by its velocity Vi. Particles collide at a constant 
collision rate F, and the colliding particles are chosen at 
random. In a collision between particles i and j, the 
velocity of each particle in the center of mass frame is 
multiplied by —a, where a is the restitution coefficient. 
The particles are coupled to a thermostat , which can be 
of two types. A stochastic thermostat [l^, [13| exerts a 
fluctuating force ^pi (t) and a velocity-dependent damping 
force —Xvi on each particle (with A > 0). During a finite 
time At > 0, short compared to F^^ and to A^^, the 
effect of the fluctuating force on a particle's velocity is 
proportional to \/Ai (see e.g. [l!]). In a Gaussian ther- 
mostat (37} , instead of applying a random force il^iit), 
A is assigned a negative value such that the term —Xvi 
acts as a velocity-dependent driving force, which rescales 
the velocities of all particles in every time interval. The 
NESS achieved with a Gaussian thermostat is exactly the 
freely cooling state, rescaled at every time step (ssj . 

For both types of thermostats, the velocity of particle 
i evolves according to 

{ value: probability: 
(1 - \At)vi + ipiVAi 1-TAt , (5) 

where Vi , Vj , and in the right-hand side are taken at 
time t. This represents the stochastic evolution of the 
system, where at probability F per unit time, particle i 
collides with a randomly chosen particle j and their ve- 
locities are updated accordingly [second line in Eq. 
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Keeping terms up to linear order in the short time in- 
terval At, the velocity of particle i is updated due to 
the interaction with the thermostat only for time inter- 
vals during which it did not collide with another particle 
in the system [first lino in Eq. ([5])]. For a stochastic 
thermostat, '4'i{t) is an uncorrelated random accelera- 
tion with (-0) = 0, and {ijP) = 2XTb, where Tb is a 
parameter, interpreted as the temperature of a bath to 
which the system is connected. For a steady-state of elas- 
tic particles {a — 1), Tq = (z;^) = Tg, as expected in 
equilibrium. For a Gaussian thermostat, ^pi(t) = and 
A = {a^ — l)r/4, which maintains a constant Ta- 
in order to obtain the velocity autocorrelation func- 
tion, we multiply Eq. fS]) by Vi{0) and average over 
stochasticity and initial conditions, to get 



v,{Q)v,{t + At)) = 

(1 - \At) {v, (0) V, (t)) + {v, (0) (t)) VAt (1 - TAt) 



TAt. (6) 



Noting that ('y^(O)V'i(t)) = 0, {v^ iO)vj{t)) = 0, neglecting 
terms of order (At)^, and rearranging, we take the limit 
At ^ to obtain 



dt 



(7) 



where kid = A + -1(1 + a)r is the "effective drag" expe- 
rienced by a particle due to the drag of the thermostat 
and the collisions with other particles. This differential 
equation has the solution 

{v{0)v{t)) = (v^)e-^^"K (8) 

Using Eq. ^ in the definition of D{lu), Eq. ([1]), we find 



D{u) 



Kid + 



(9) 



To calculate the mobility, a periodic acceleration ^(t) = 
^oe''^' is applied to a single particle i. Due to the mean- 
field assumption, and in the thermodynamic limit of a 
very large system, the probability distribution of Vj(t) 
{i j) is unaffected by the force ^{t) acting on particle 
i. From Eq. ([5]), we similarly have 



dt 



= -Kin {vi{t)) +^oe 



(10) 



This has the steady solution {vi{t)) = ^oe*'^V(^ + 
from Eq. ([2]) we see that 



1 



Kid + 



(11) 



Using Eqs. (H]), ©, and (fTT|) we find that the equilibrium 
GK formula holds with Tpd{uj) = Tq, even though the 
system is far from equilibrium. 



We now turn to the derivation of the GK relations in 
the Maxwell model in higher dimensions. We will dis- 
cuss the three dimensional case in what follows, but the 
derivation applies to any dimension, and generalizes pre- 
vious ID [20, 22,] and 2D ^ results. 

The equation of motion for the velocity in the lab- 
oratory frame is 



value: 



V. (t + At) 



(1 - XAt) V, - 
V, - i±2 (n • 



Av 



At 
y)n 



probability: 
1-TAt 
TAt 



(12) 



where the quantities in the right-hand side are taken at 
time t. As in the ID derivation, A denotes the drag co- 
efficient, the vector xpi is the random driving accelera- 
tion applied to particle i, and a is the restitution coef- 
ficient, n = (ri — rj) /2R is the unit vector between the 
centers of the colliding particles. R is the particle radius. 
Avy = Vj — Vj is the precollisional relative velocity. 

As in the ID case, the derivation of the GK relation 
consists of separate calculations of the mobility and au- 
tocorrelation; the GK relation is then obtained by taking 
the ratio of the two results. 



A. Mobility 

To calculate the mobility, we add an oscillatory driving 
to Eq. dm): 



V, {t + At) 



value: 



probability: 



4oe"^*Ai+ { (1 - XAt) V, + i(>^y/Ai 1 - TAt 
V, -i±ii(n-Av,y)n TAt 



(13) 



Averaging over the stochasticity and the initial distribu- 
tion we find: 

(v, (t + At)) = Coe'"*At +{1- TAt) (1 - AAt) (v, (t)) 

+ TAt (v, (t)) - TAti^ ((fi • Av„ ) n) , 

(14) 

where the term involving ipi vanishes, since (ipi) ~ 0. 

Consider the collision of two particles, i and j. Fig- 
ure [1] depicts this collision, in the reference frame moving 
with the pre-coUision velocity of particle i. b is the di- 
mensionless impact parameter, < 6 < 1, distributed 
with probability (c? — 1)5'^^^ in d dimensions, for d > 1 
(for d = 1, b = identically). Due to the assump- 
tion of molecular chaos, 6 is independent of the prec- 
oUisional velocities Vj (0) , v, {t) ,Vj (t) and Av^^ (t). Let 
us decompose n = n,, -I- n_L, such that n,, is aligned with 
the colhsion direction Avy {t), and perpendicular to 
that direction, see Fig. [TJ Note that n,!,!!^ are not of 
unit length. Now define k||,kj^ to be the unit vectors 
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B. Autocorrelation 
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FIG. 1: Collision of two particles in the rest frame of particle 
i before the collision. 



in the nii,n^ directions respectively. Then nj 



6k I 



and n,! = \J\ — 6^k||. Substituting these definitions and 
neglecting terms of order (At) , the preceding equation 
reads 

(v, (< + At)) = ^oe'"* At + (1 - AAt) (v, (t)) 
1 + a 



rAti±^((l-6^) (k„.Av,)k, 
rAtl±^((n„.Av,,)n^). (15) 



Since k,, and Avy are co-aligned, 



k, • Av,;, k„ = A 



and because h is independent of Avy (t), 



\-h'\ k,- A 



(16) 



= ((l-62)Av,,.) = (l-62)(Av,,). (17) 

The term ((n,, • Avij)n^) in Eq. ([15]) vanishes (without 
assuming it can be decomposed into uncorrelated terms), 
since it is anti-symmetric in n_L: for every n^, the contri- 
bution of (n,i • Avy ) to the average, with n_L — > — nj^, 
is equally probable, and with opposite sign. Hence, using 
Eq. ((ni), Eq. becomes 

(v, (t + At)) = ^oe'"*dt + (1 - AAt) (v, (t)> 

-rAti±^(l-52)(Av,,) 

= 6e'"*At + (1 - AAt) (v, (t)) 
1 + ,2\ 



TAt- 



(l-62)(v,(t)), (18) 



since (v^ (t)) = 0. Gathering the terms, we obtain the 
differential equation 



dt 



= -k{v, (t))+|oe^"*. 



with n = \ + ^{l-P)T. Hence 

1 



K + iuj 



(19) 



(20) 



as in the ID case, but with a dimension-dependent effec- 
tive drag coefficient k. 



Using Eq. ([12]), we compute the change in (v^ (0) (t)) 
during a short interval At. We multiply both sides of 
Eq. (fT^ by (0), and take the ensemble average. Noting 
that {ipi (t) v; (0)) = 0, and n_L • Avy (t) = 0, we find 

(v, (0) V, (t + At)) = (1 - AAt) (v, (0) V, (t)) 



EAt- 



((n„.Av,,)[(n„ + n^)-v,(0)]). (21) 



With k,, as defined abc 



this becomes 



(v, (0) V, (t + At)) = (1 - AAt) (v, (0) V, (t)) 
l + a 



-FAt- 



1-b' 



k, • A^ 



^^ (0) 



-FAt- 



1 



((n„ 



Av,,)[ni -v, (0)]) 



(22) 



2 

As in the mobility calculation, the term 
((n,! • Avij) [n_L • Vi (0)]) vanishes, because it is anti- 
symmetric under the transformation: — nj^. 

Moreover, ^k„ • Av^^ k,, = Av^ . Rearranging we find 



d(v, (0).v, (t)) 
dt 



with K as defined following Eq 
ID case, 

D{lo) = 



K (v, (t) • V, (0)) (23) 
T9l). Hence, as in the 



J- (24) 

Eq. ([M]) and ([^ are generalizations of Eq. ^ and 
([TT]) respectively; Indeed, in ID, 6 = 0, and the result 
i^iD = \ + \T {1 -\- a) is reproduced. In 2D and 3D, 
K2D = A + |F (1 + q), and n^o = A + |F (1 + a). As in 
the ID case, Eq. ([20]) and ([24]) show that the mobility 
and autocorrelation functions have the same frequency- 
dependence, hence the frequency-dependent GK relation 
holds for the Maxwell model, with Tpd(lo) = Tq, in any 
dimension, and for both customarily studied thermostats. 



IV. MEAN-FIELD GRANULAR GAS 

In this section the GK relation is discussed for the 
MFGG model which has a velocity dependent collision 
rate. The zero frequency GK relation (Einstein relation) 
for the 3D version of this model has been studied previ- 
ously. Using kinetic theory approximations, the value of 
TpD (0) was shown to be close to, but not exactly equal 
to Tg- For a tracer particle of the same mass and size as 
the rest of the particles, Tpo (0) deviates from Tq by up 
to 1% for the stochastic thermostat and by up to 6% for 
the Gaussian thermostat [l3| . 

We now show that in the infinite frequency limit the 
FD temperature is exactly equal to the granular temper- 
ature: 



Tg 



1. 



(25) 
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In fact, what we show is that for frequencies much higher 
than the highest characteristic frequency in the system, 
this relation holds. Eq. ([25]), together with the results of 
[13], show that the FD relation is violated. They do, how- 
ever, indicate that the violation, measured as the varia- 
tion in Tfd (i^) with frequency, is small. 

Consider first the mobility The mobility in the 
Maxwell model was calculated using Eq. . This equa- 
tion must be modified for the MFGG. This is because the 
effective drag experienced by a particle depends on its in- 
stantaneous velocity: k, = k{v). Note that, due to the 
mean-field (molecular chaos) assumption, the collisions 
of a particle do not depend on the history of the parti- 
cle, and therefore k only depends on the current particle 
velocity. The analog of Eq. (|19p can therefore be written 
as 



d{^^ it)) 

dt 



= -{k (v,) V, {t)) + i (t) 



(26) 



where as before the function ^ (t) is the force applied to 
the tracer particle. We denote the velocity distribution 
of the tracer particle at time t hy p (v^ , t) . Writing the 
averages explicitly, Eq. ([26]) becomes 

J -^^^^^Vidv^ ^ ~ J p{vi,t) K{vi) Vidwi + ^ (t) 



and its Fourier transform reads 

K (Vi) 



iuj / p(vi,w) 



1 



(27) 

v.rfv, = , (28) 



where p (vi, w) , ^ (w) are the Fourier transforms of 
p (vi, t) , ^ {t) respectively. Rearranging, and noting that 
(vi [ld)) and ^ (w) are vectors in the same direction, so 
that their quotient is well-defined, we find that 



Pl[uj) 



(v^ (^)) _ / pj^i.Uj) VtdVi 



o 



^(29) 

where O (^) denotes a function of order {l/i^j) , for 
l/w^O. 

To calculate the fluctuation D{uj), we denote: 
C (t) = (v(O)v(t)), and the Fourier-Laplace transform 
of a function by FL [. .] . Then: 



FL [C{t)] 





-1 
ito 
1 



C{t)e- 



C{t)e- 



*dt 



it=o 



dC jt) 
dt 



C (0) + FL 



dC {t) 



dt 



'dt 



(30) 



This recursion relation [39| is satisfied by the infinite se- 
ries 



n— ^ ' 



dt"- 



t=o+ 



FD 



Tg 



Maxwell model 
mean-field granular gas 
schematic interpolation 



UJ 



FIG. 2: (Color online) FD temperature in the Maxwell model, 
and in the mean-field granular gas model. In the former, the 
temperature obtained from the Green-Kubo relations is fre- 
quency independent and equal to the granular temperature. 
In the latter, the temperature obtained depends on the mea- 
surement frequency. The schematic curve interpolates be- 
tween the high and low u limits. The difference between the 
two curves is not to scale. 



If all the derivatives of C (t) are well-behaved, we have 
D{cj)^~C{0) + o(^]. (32) 



Using Eq. ^ and 
is 



the resulting FD temperature 



^-ic(0)4(v^)^T„ 



(33) 



(31) 



which is what we wanted to show; c.f. Eq. ([25)) . 

This result can be interpreted as follows: the system 
possesses two time scales, the mean time between colli- 
sions, and the thermostat time scale, A~^. For frequencies 
higher than the inverse of these two time scales, D (ui) 
and fj, (lu) have the same frequency dependence, and the 
FD relation holds, with D (uj) /Tg ~ /i (ut) ~ ■ 



V. DISCUSSION 

The results for the two models can be summarized as 
follows, see Fig. [21 In the Maxwell model the FD relation 
holds, with a single FD temperature for all frequencies. 
In the more physicall y re alistic MFGG, the FD relation 
is violated: whereas gives Tpu (0) ^ Tq, we have 
found Hhi^^qoTfo (uj) = Tg, demonstrating that Tpu 
is not independent of lu. The difference in Tpo between 
the limits of high and low frequency is smaller than sev- 
eral percent. We would thus expect that the variation 
in Tpo (w) is small for all frequencies. This variation is 
expected to be the largest for small a and could perhaps 
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be detected in simulations. Molecular dynamics simu- 
lations are hard to conduct for very inelast ic p articles 
(a w 0) , since the grains tend to cluster [iO, El] . In 
direct simulation Monte Carlo techniques, on the other 
hand, the mean-field property is explicitly imposed, and 
the variation of Tpo with frequency may be easier to 
detect. It would furthermore be interesting to explore 
whether analytic progress may be achieved in studying 
the high frequency limit for dense systems where corre- 
lations give rise to significant FD violations [l^]. Last, 
we note that the technique we used to calculate the in- 
finite frequency limit for the MFGG model applies to a 
more general class of models, such as models in which 



the restitution coefficient depends on the particles' rela- 
tive velocity [ii,!!!!!!. 
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